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SERIE |. CIENCIAS, NO 38 Pe ISSN 0037-8437 


FOTOLISIS UV DE 5-BROMURACILO EN SOLUCION ACUOSA 
INFLUENCIA DE OXIGENO Y D-—DESOXIRIBOSA (1) 


Por E. GILBERT y C. CRISTALLINI E 


UV-Photolysis (254 mm) of 5-—bromouracil in aqueous solution in the presence of oxygen was 
investigated at pH 6. Isodialuric acid and 4,5-dihydroxy-imidazolidone-2 as a secondary photuproduct 
were identified. The bromine atom from the homolytic splitting of the C-Br-bond and molecular bromine 
have no effect on the product distribution under our experimental conditions. In the presence of 
2-deoxy-D-ribose and oxygen the UV-photolysis of bromouracil leads to uracil, isodialuric acid and 
4 5-dihydroxy imidazolidone-2. The yield of uracil increases with increasing deoxy-D-ribose concentration 
while the yield of isodialuric acid and 4,5-dihydroxy-imidazolidone decreases. 


La introducción del S-bromouracilo (Bu) en lugar de timina en el ADN de fagos y 
bacterias aumenta su sensibilidad a la radiación ultravioleta.* Como contribución para 
interpretar el mecanismo de acción, resulta de interés investigar la influencia de oxígeno y de 
2-d-desoxiribosa sobre la fotólisis-UV del Bu. 


RESULTADOS Y DISCUSION 
a. Fotólisis-UV de bromouracilo en presencia de oxigeno 


La fotólisis-UV(A = 254 nm) del Bu (c = 5.10% mol/1) en solución acuosa, saturada con 
oxígeno (pH = 6, sin buffer, 0, = 1 atm) conduce fundamentalmente a dos productos 
orgánicos de fotólisis: el ácido isodialúrico (Id) y la 4,5-dihidroxiimidazolidona-2 (Im) (Fig. 1 
E 

Como se observa en la ¡igura 2, inicialmente se forma sólo 1d; manteniendo la exposi- 
ción el Id decrece y se forma como producto secundario Im. Con referencia al mecanismo de 
formación de Id se supone que la etapa fotoquímica inicial es la ruptura homolítica de la 
ligadura C—Br (1). 

La reacción del radical uracilo con el oxígeno para dar Id ha sido previamente estudia- 
da? en la fotólisis del iodouracilo (Tu). 

En el presente trabajo se admiten mecanismos de formación idénticos para Id. En la 
fotólisis-UV de lu se encontró hasta el presente solamente 1d?. Debido a los rendimientos 
cuánticos mucho menores de Bu en comparación con lu (véase parte experimental), la 
fotólisis del Bu requirió tiempos de exposición diez veces mayores y concomitantemente una 
reducción considerable de la distancia entrela fuente y la solución para que la conversión 
fuera mensurable. Usando tiempos de irradiación prolongados se encontró que el Id formado 
inicialmente. se convierte luego fotolíticamente en la Im. Aún agregando Id (c = 5.10% 
mol/1) previamente a la irradiación a la solución de Bu, este no pudo ser detectado a los 90 
minutos de exposición por cromatografía en capa delgada. Se demostró también que con una 
exposición mayor de 90 minutos de lu, el Id formado durante la irradiación se convierte en 
Im por un proceso secundario (Fig. 3). 


* Traducción del trabajo publicado en Z. Naturforsch. 28 b, 615-619 (1973): “UV-Photolyse von 5-Bromou- 


racil in wássriger Lóusung”. 


2 Instituto Nacional de Tecnología Industrial. 
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La fotólisis de Id no ha sido estudiada exhaustivamente. Hasta el momento hemos 
obtenido los siguientes resultados: A partir de experiencias con Bu.marcado con *C en 
posición 2 y 4 respectivamente, hemos podido establecer que el átomo de C en posición 4 se 
desprende. El radiocromatograma de la solución irradiada de Bu4-'* C (Fig. 4) muestra como 
productos activos únicamente ld y el Bu sin reaccionar, en cambio no se observa la 
formación de Im como en el radiocromatograma de la solución irradiada de Bu-2-** C (Fig. 1) 
Además en la solución irradiada de Bu-2-** C se detectó '* CO, lo cual significa que durante 
la descomposición fotolítica el Id sufre decarbonilación dando un anillo de 5 átomos (Im). 
Wake y colaboradores observaron también la formación de un anillo de 5 átomos en la 
fotólisis-UV del aloxano?. Actualmente no podemos dar más detalles acerca del mecanismo 
de formación de Im. ] 

Queda aún por discutir la reacción del átomo de bromo separado de acuerdo con (1). 
Por dimerización podría formarse bromo molecular, pero este no fué detectado en la 
solución irradiada. Además, como muestran los ensayos en blanco (sin irradiación), el 
bromouracilo reacciona cuantitativamente con bromo formando 5,6-dibromouracilo (Dbu), 
(1). 

El Dbu es estable y puede ser aislado de la solución acuosa (véase parte experimental). 
Sometido a la fotólisis se forman 4 productos no identificados. Ni el Dbu ni sus productos 
de fotólisis son detectados en la fotólisis UV de Bu. 
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En otros ensayos realizados se hizo gotear agua de bromo (1,5 ml, e = 2.10* mol/1) a la 
solución del Bu (c = 1,5.10-3 mol/1) durante la irradiación de una hora. Tampoco en este caso p 
pudo ser detectado bromo molecular. Como productos de fotólisis resultaron únicamente 
Id y Im; en la misma relación se formaron Id y Im en un ensayo en el cual se dejó gotear 
agua durante la irradiación. En un ensayo paralelo se hizo gotear agua de bromo sobre agua 
pura que fué irradiada. Después de la irradiación no pudo ser detectado bromo molecular. En 
base a los resultados obtenidos se concluye que el bromo separado homolíticamente no se 
adiciona al Bu. Se supone que el bromo se oxida durante la irradiación y que en su mayor 
estado de oxidación no influye en la fotólisis de Bu ni sobre la formación de los productos 
de fotólisis. 
Existe aún la posibilidad de eliminar el bromo después de la homólisis de la ligadura C-Br. 
Esto se logró mediante IK, que se agregó a la solución de Bu (eo L0% mol/1) en una 
concentración de 5.10% mol/1. 

En la fotólisis de Bu en presencia de ioduro se forma ¡odo libre, resultante por una parte 
de la reacción de radicales peroxiuracilos con ioduro * (IV). 
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Por otra parte según (1), los átomos de bromo separados Ó el bromo molecular resultante de 
su dimerización reaccionan con el ioduro (V), (VI). 


A Bro da (V) 
Br) > BD L (VI) 


Admitiendo que las reacciones 4 y 5 Ó 6 son cuantitativas se formaría un mol de iodo por 
cada molécula de Bu consumida. Esto se cumple únicamente con tiempos de irradiación 
cortos (Fig. 5). En el transcurso subsiguiente de irradiación la formación de iodo es menor 
que lo que correspondería a la formación de Bu. Esto puede atribuirse en ciertas condiciones 
al hecho que la concentración de ioduro disminuye y que las reacciones 4 y 5 Ó 6 pueden no 
ser cuantitativas. 


b. Fotólisis-UV de Bromouracilo en presencia de oxígeno y 2-d-desoxiribosa. 


En la fotólisis-UV de bromouracilo en solución acuosa en presencia de oxigeno y 
desoxiribosa (DR) se encontró según la concentración inicial de la DR los siguientes 
productos en relaciones variables: uracilo (U), Id e Im. La cantidad transformada de Bu 
corresponde a la suma de los productos de fotólisis. Dado que en ausencia de DR no se 
forma uracilo, se admite que el radical uracilo formado según (I), en presencia de donores de 
hidrógeno, simultáneamente con la adición de oxígeno (II), da lugar a otra reacción que 
conduce a la formación de uracilo. Esto sólo puede efectuarse mediante una abstracción de 
H en la DR(VII) 
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Reacciones del radical U con otros donores ya han sido descriptas previamente*. Si las 
reacciones 2 y 7 fueran competitivas, la relación de los productos U a Id e Im resultarían 
desplazadas a uracilo al disminuir la concentración de oxígeno manteniendo la concentración 
de la DR constante. Esto fué corroborado por fotólisis de soluciones de Bu a saturación de 
oxigeno y aire a concentraciones constantes de DR (Fig. 6). 

Estos resultados permiten afirmar que la etapa fotolítica inicial de la fotólisis-UV de Bu 
consiste en una ruptura homolítica de la ligadura C-Br, como ya fué formulado por otros 
autores”. 


PARTE EXPERIMENTAL 


Las condiciones experimentales no explicitadas son las mismas que en el trabajo prece- 
dente”. 
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Condiciones de irradiación 


Como fuente de UV se usó una lámpara de Hg de baja presión (A. Grántzel, Karlsruhe). 
El flujo cuántico a 254 nm (4,9.10*% cuantos/min) fue medido según método de Hatchard y 
Parkeró. Se elimina la radiaciónde 185 nm por filtración con placa de vidrio Vicor. Las 
soluciones a irradiar fueron colocadas en celdas de 1 cm de paso óptico en un block 
termostatizable directamente delante dela lámpara. En los ensayos con lu se irradió a una 
distancia de 2 cm de la lámpara. Los rendimientos cuánticos se calculan para Bu en: 1,5.1073 
y para lu en: 2,1.10-2. 


Bromo 


El bromo se detectó mediante reacción con ioduro. Se utilizó como solución de ioduro 
una solución buffer según Hochanadel”. 


Uracilo, ácido isodialúrico 


La determinación cuantitativa de uracilo y Id se realizó por radiocromatografía. Se 
identificó por comparación de los valores Rf con productos sintéticos. Los datos cromatográ- 
ficos de Id y U han sido descriptos previamente?-8. 


Determinación de CO 


El gas formado en la solución irradiada de Bu-4-'* C se barrió con corriente de N, 
haciéndolo burbujear en una solución de Ba(OH),. Anteponiendo a la solución de barita una 
solución de PdCl,, que oxida el CO, puede determinarse la actividad. Sin interceptar con 
solución de PdCl, la actividad medida resulta 5 veces menor. 


4,5 Dihidroxiimidazolidona-2 


Se aisló Im a partir de 2 litros de una solución irradiada de Id(c = 10? mol/1). En un 
dedo de cuarzo se irradió cada vez 70 mi de solución de Id mediante 2 lámparas de Hg de 
baja presión durante 90 minutos. La solución irradiada se concentra a 300C. A partir de la 
solución concentrada cristalizó Im al cabo de 14 días a OC. Con fines comparativos se 
sintetizó Im a partir de urea y glioxal según síntesis de Reibnitz?. El Im preparado por 
fotólisis (Im fot) y el sintetizado (Im sint) responden al mismo espectro de absorción en el 
infrarrojo. Eos espectros de resonancia magnética nuclear de *H, en (CD, ),SO, del Im fot. y 
del Im sint. coinciden. Ellos muestran como es de esperar tres señales de la misma 
intensidad: 7= 5,86 (protones de OH, doblete, J = 6,0 Hz) y 7= 4,63 (protones de CH, doblete, 
J= 6,0 Hz). 

El peso molecular determinado por osmómetro de presión de vapor fué para el Im fot. 
de 114 + 4, siendo el valor teórico de 118. 
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Analisis de Im fot. 
Calculado C 30,000 HS5,l N 23,8 O 40,51 
Hallado C13068. H:52 N 23,67 040,1 


Im fot y Im sint muestran los mismos valores de Rf en distintos medios eluyentes (Tabla 1). 


TABLA 1: Valores de Rf obtenidos sobre capa delgada de celulosa para el Im obtenido por 
fotólisis de Id y para el Im sintetizado en este laboratorio. 


No Eluyente Valor Rf Valor Rf 
Im fot. Im sint. 
1 Acetato etilo/ác. fórmico/H, 0;10/2/3 0,28 0,28 
2 t—butanol/H,O/ác. acético; 44/11/0,26 0,47 0,47 
3 Acetato de etilo/ác. fórm/H,0O; 50/20/20 0,44 0,44 
4 isopropanol/H) 0; 7/3 0,57 0,57 


La detección sobre la capa cromatográfica fué realizada por vaporización con solución 0,1 N 
de NaOH. Después de 20 minutos se forman manchas verde-oscuras sobre fondo blanco. 


5,6 Dibromouracilo 


Se prepara 5,6 dibromouracilo por la reacción de bromo con bromouracilo. A una 
solución acuosa de Bu 10% molar se dejó gotear bromo a temperatura ambiente en relación 
1;1. La solución queda incolora, es decir que el bromo reacciona cuantitativamente. El valor 
de pH disminuye de S a 3. El bromuro se determinó cuantitativamente con NO3Ag. 
Concentrada y enfriada la solución a O”C se obtuvo un producto cristalizado. El análisis del 


producto bromado y el peso molecular corresponden a un dihidrato de 5 6-dibromouracilo 
(Dbu). 


Andlisis de Dbu.2 H,0 


Calculado EMULST H 1,96 N 9,14 O 20,9 Br 52,2 
Hallado E-16:05=.H:1:91 N 8,1 O 1964 Br 53,9 


Determinación de peso molecular: 305-5 (teór. 306) 
en la Tabla 11 se muestra los valores Rf con distintos eluyentes. 


TABLA || 


VALORES DE RF PARA 5,6-DIBROMOURACILO OBTENIDOS 
SOBRE CAPA DELGADA DE CELULOSA 


Eluyente Valor Rf 
1 0,83 
Z 0,86 
3 0,87 
4 0,81 
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-2 


LMP.1O0 


—= 


0,2 0,4 0/6 : 0,58 


—> varon De Rf 


Fig. 1. — Radiocromatograma de una solución irradiada de S-bromó uracilo-2-* C, revelado sobre capa 


delgada de celulosa con eluyente II (Tabla 1). tiempo de exposición: 60 minutos (sobre condiciones de 
exposición ver varte experimental). 


E 10% (mo111) 


————>CONC 
6 
o 


20 40 60 80 100 120 


———> TIEMPO DE IRRADIACION t (min) 


Fig. 2. — Fotolisis-UV (A = 254 nm) del S-bromouracilo (C = 5.10% mol/1) en solución acuosa (pH = 6,0, 
x l atm) U =Bu transformado, €8 = ácido isodialúrico, xx = 4,5-dihidroxi-imidazolidona-2. 
á 


CONC. 1 0 (me111) 


20 40 60 80 100 


———>> TIEMPO DE IRRADIACION t(min) 


Fig. 3. — FotólisisUV (A = 254 nm) del lu (c = 5.10% mol/1) en solución acuosa (pH = 6,02 = 1 atm.). 
69 — ácido isodialúrico, xx: 4,5-dihidroxiimidazolidona-2. 
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Fig. 4. — Radiocromatograma de una solución irradiada de S-bromouracilo-4-14C 
mentales son análogas a las dadas en la Fig. 1. 


- Las condiciones experi- 
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A Tiempo de irradiación t (min) 


Fig. 5. — FotólisisUV (A = 254 nm) del Bu(c = 5.104 mol/1) en solución acuosa en presencia de IK (c = 
5.10%mol/ 1) OU = Bu transformado por fotólisis; € = Id; O: Im; Á=1L. 


—— conc. 10*(mo//1) 
o 
u 


-3 -2 a 


——>lo0g C 2-Desoxi-d-ribosa (mol/1) 


Fig. 6. — Fotólisis-UV de Bu(c-5.107? mol/1) en presencia de DR y oxígeno 


: solución saturada de 
oxígeno; ——: solución saturada de aire; O) = uracilo; 6 = 2 Idy Im. 


Nota: los valores dados de los distintos productos han sido tomados a partir de diagramas de depedencia 
de dosis cuando ha reaccionado el 20% de Bu. 
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BORRADO DEL EFECTO DE MEMORIA EN MEZCLAS DE CRISTALES 
LIQUIDOS NEMATICOS Y COLESTERICOS 


Por BETTY KERLLEÑEVICH! y ANDRE COCHE? 


RESUMEN 


Se analizaron las características del proceso de borrado por aplicación de una tensión alternada de 
audiofrecuencia, en mezclas de distintas proporciones de cristales líquidos nemáticos y colestéricos. Se 
determinó la influencia de varios parámetros sobre el tiempo de borrado: 1) temperatura (y viscosidad), 
2) valor de la tensión alternada de borrado, 3) concentración de producto colestérico, 4) campo continuo 
que produce la difusión dinámica y efecto de introducir un retardo entre su supresión y la aplicación de la 
tensión alternada. 


INTRODUCCION 


Los efectos Ópticos de memoria en mezclas de sustancias nemáticas y colestéricas fueron 
descriptas por primera vez por Heilmeier y Goldmacher1.2. Estos autores demostraron que 
una mezcla de 90% en peso de anisilidén-para-aminofenilacetato (APAPA) y 10% de nona- 
noato de colesterilo (NC) produce una intensa difusión de la luz, en forma análoga al 
fenómeno de difusión dinámica (dynamic scattering) observado en los productos nemáticos 
puros, cuando se la somete a un campo eléctrico continuo o alterno de baja frecuencia 
(< 50 Hz). Pero en el caso de las mezclas, el efecto persiste después de suprimido el campo. 

El “borrado” de la información así almacenada puede lograrse aplicando a la célula una 
tensión de audiofrecuencia superior a 500 Hz. Los mismos autores determinaron que el 
tiempo de borrado varía linealmente con V-1 donde V es la tensión de borrado. Varios 
autores3.4 demostraron que la aplicación de campos eléctricos de frecuencias altas da como 
resultado la formación de la estructura colestérica de planos de Grandjean, cuya claridad 
Óptica es muy parecida a la del producto nemático, mientras que los campos de baja 
frecuencia producen una estructura focal-cónica: su apariencia “policristalina” produce una 
difusión intensa de la luz, que es más notable cuanto más gruesa es la película estudiada. 

Vamos a presentar algunas características del proceso de borrado de información almace- 
nada por aplicación de un campo eléctrico continuo, en mezclas de LCI (producto nemático 
para el rango de temperaturas entre 5 y 70%C)? y de nonanoato de colesterilo (NC). 


MONTAJE EXPERIMENTAL 


Las células que 'hemos utilizado en nuestras experiencias? constan de una película delgada 
de la mezcla nemática-colestérica entre dos electrodos transparentes, en nuestro caso de 
vidrio revestido de una capa de óxido de estaño, entre las que se disponen separadores de 
mylar. La distancia entre electrodos se mantiene uniforme mediante un dispositivo de presión 
constante. Para obtener un ordenamiento lineal y homogéneo de la mezcla de cristales 
líquidos, se frotan previamente los electrodos transparentes sobre papel siguiendo la técnica 
de Chatelain. Los electrodos pueden estar a distancias que hemos elegido de 18 um y en 


1 Departamento de Ciencias Exactas, Universidad Nacional del Sur, Bahía Blanca. 
2 Centre de Recherches Nucléaires, Strasbourg, Francia. 
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algunos casos de 36 um, y éste es en consecuencia, el espesor de la película de cristal 
líquido. La célula se ubica en una caja termostatizada a + 0,1%C en la que la temperatura 
puede hacerse variar entre 10 y 50*C aproximadamente. 

La célula se ilumina con luz blanca y la luz que transmite se hace incidir en un 
fotomultiplicador; la señal anódica de éste se introduce a un osciloscopio o a un registrador 
X-Y. De este modo puede determinarse a cada instante la intensidad de la luz transmitida I, 
y la de la luz difundida I,=I,—I,, siendo I7 la intensidad transmitida por la célula en 
ausencia de campo eléctrico. De las curvas de decrecimiento de la luz difundida, puede 
deducirse el tiempo de borrado 7, que definimos como el tiempo en que la luz difundida Is 
decae de 90 a 10% de su valor inicial, cuando simultáneamente con la supresión del campo 
continuo se aplica un campo alterno. Hemos adoptado como frecuencia de la tensión de 
borrado 1 KHz ya que de experiencias preliminares se comprobó que el tiempo de borrado 
varía muy poco con la frecuencia, entre 0,5 y 5 KHz. 

Hemos estudiado las variaciones de Tp con varios parámetros: temperatura, tensión de 
borrado, tensión continua aplicada, para diferentes concentraciones de NC que se variaron 
entre 1 y 9% en peso. Además se trazaron las características de corriente en función de la 
tensión, para las distintas temperaturas, en cada célula analizada. 


RESULTADOS 


a) El agregado de sustancias colestéricas a las nemáticas, aún en pequeñas proporciones, 
produce la aparición de ciertas características de las sustancias colestéricas y en particular 
el efecto de memoriaf. En realidad esos efectos aparecen a partir de una tensión V y 
superior a la tensión umbral V,, de la difusión dinámica para la sustancia nemática pura. La 
tensión V,, corresponde aproximadamente al punto de discontinuidad de la curva de 
corriente en función de la tensión de la mezcla. En un trabajo anterior” (Fig. 1) hemos 
observado que su valor aumenta con la concentración de NC y disminuye para temperaturas 
más elevadas. Si la tensión aplicada tiene un valor comprendido entre V,, y V,,, la difusión 
dinámica desaparece espontáneamente al suprimirse el campo continuo, según una constante 
de tiempo cuyo valor no difiere apreciablemente de la del producto nemático en las mismas 
condiciones. 

b) Para un espesor dado de la célula, al hacer aumentar la temperatura dentro del rango 
de 10 a 50*C, el tiempo ”p decrece. En Fig. 2 se ve un gráfico semilogarítmico de esta 
variación. Fuera de los valores de Tp obtenidos para temperaturas del orden de 40 a 50*C, 
el decrecimiento de Tp sigue una ley similar a la observada para el decrecimiento espon- 
táneo en el caso del producto nemático puro. Esta ley es del tipo 


0EÉ 
kT 


(1) 


TD = Á exp 


donde k es la: constante de Boltzmann y A es una constante arbitraria. La cantidad 6E=0,4 
eV (determinada a partir de las curvas) no varía apreciablemente con el aumento de 
concentración de NC entre 1% y 9% en peso. Este valor es del mismo orden que el 
determinado para el producto nemático LCI puro. Para estas concentraciones la variación de 
Tp con la temperatura parece relacionada directamente con la viscosidad n del producto 
nemático. En los trabajos realizados con LCI puro, llegamos a determinar que la variación de a 
con la temperatura se debe exclusivamente a la de la viscosidad de la sustancia nemática, y estos 
resultados coinciden con los obtenidos por otros autores, Koelmans y Van Boxteló. quienes 
estudiaron el tiempo t, de crecimiento de 1, en mezclas de bases de Schiff. La variación 
de Tp (y de n) sigue prácticamente la.misma ley, dada por ec. (1) Ver Fig 3. 
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El viscosímetro utilizado es similar al descripto por Lindermann y Rikmenspoel? y 
permite efectuar mediciones sobre pequeños volúmenes (<1 cm? ). El principio se basa en el 
ocultamiento sucesivo de dos orificios, por una bolilla de acero que cae dentro de la 
sustancia cuya viscosidad es estudiada. La luz que normalmente pasa a través de los agujeros 
para incidir en un fotomultiplicador es cortada entonces por el paso de la bolilla. Resultan 
así dos pulsos sobre la señal de salida observada en el osciloscopio o sobre un registrador 
X-Y. El tiempo que separa los dos pulsos depende de la viscosidad. 

c) Nuestros resultados experimentales muestran que el tiempo 7, varía con la tensión 
alterna de borrado V, siguiendo la relación: 


donde a es una constante. Se han representado ejemplos de esta variación en Fig. 4 en escala 
doble logarítmica, para varias concentraciones. Se verifica dentro del rango de concentracio- 
nes estudiado que el exponente n es prácticamente independiente de la temperatura y del 
campo continuo que produce la difusión de la luz; crece siguiendo una ley lineal con la 
concentración de NC como se ve en Fig. 5 y su valor es, por ejemplo: n= 1,5 y 3,15 para 
1% y 8,5% de NC, respectivamente. 

d) La Fig. 6 muestra que para una temperatura dada y para una misma tensión de 
borrado, el tiempo 7, crece linealmente con la concentración de NC entre 1% y 5 (0 6)% 
de NC. Para concentraciones mayores, este tiempo de borrado aumenta mucho más rápida- 
mente. 

e)Se estudió también la influencia sobre el tiempo de borrado Tp del campo que 
produce la difusión de la luz. Para concentraciones de NC superiores al 4%, el tiempo 7, 
decrece considerablemente cuando la tensión continua se hace superior a los 20 Volt (en el 
caso de una célula de 18 um de espesor) y cuando la tensión de borrado se aplica en el 
mismo instante en que se corta el campo continuo. Esta variación se observa en Fig. 7. En 
cambio para concentraciones menores Ty es prácticamente independiente de la tensión 
continua. Este efecto podría deberse a cambios de textura en el cristal líquido inducidos por 
el campo. Se trataría de transiciones de estructura focal-cónica a Grandjean, que pueden 
explicarse mediante consideraciones energéticas y que reducirían el efecto de memoria y el 
liempo Tp :"> 

Hemos observado también que si la tensión de borrado se aplica al cabo de un tiempo 
Tp después de suprimido el campo continuo, 7, aumenta en consecuencia, y para un 
tiempo Tp suficiente (dos minutos aproximadamente), tiende hacia el orden de valores 
correspondiente a las bajas tensiones continuas (20 Volt para el caso de la célula de 18um). 
Este efecto podría explicarse por una transformación de textura con transición por una 
estructura más estable, cuyo mecanismo está siendo estudiado. Hemos observado el mismo 
efecto para concentraciones mayores de NC. 


CONCLUSIONES 


La degradación de la estructura almacenada por aplicación de un campo continuo a una 
mezcla de cristales líquidos nemático y colestérico, depende de varios parámetros: tempera- 
tura, tensión de borrado, concentración de producto colestérico, campo continuo aplicado. 

El tiempo que caracteriza esta degradación se hace más corto con el aumento de 
temperatura, tensiones de borrado y campo continuo. 

Para concentraciones pequeñas de producto colestérico (=1%), se observa un efecto de 
memoria y en estos casos pueden obtenerse tiempos de borrado cortos. 
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Fig. 1. — Curvas características de corriente en función de la tensión a— para algunas concentraciones de 
NC en LCI: b— para algunas temperaturas. 


Fig. 2. — Tiempo de borrado Tp en función de la temperatura para tres concentraciones de NC en las 
siguientes condiciones: tensión continua 50 Volt, tensión de borrado 80 Volt (en valor eficaz). 


Fig. 3. — Variación del tiempo de borrado Tp y de la viscosidad N en función de la temperatura, para el 
LCI puro, y variación de N para una mezcla de LCI+ 6% de NC. 


Fig. 4. — Tiempo de borrado 7;¡) en función de la tensión alterna, para varias concentraciones de NC en 
las siguientes condiciones: temperatura 240C, tensión continua 50 Volt. 


Fig. 5.— n en función de la concentración de NC en LCI (correspondiente a la ley Tp =a V *) 


Fig. 6. — Tiempo de borrado Tp en función de la concentración de NC para una tensión continua de 50 
Volt. 


Fig. 7. — Tiempo de borrado en función de la tensión continua para una concentración de 9% de NC en 
LCI y a la temperatura de 230C. 


Fig. 8. — Tiempo de borrado Tp en función del retardo Tp entre la supresión del campo continuo y la 
aplicación de la tensión alterna de borrado, para varias tensiones continuas. Concentración: 8% de NC 
en LCI, temperatura 230C. 
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“EXPERIMENTACION CON NUTRIENTES: LA DOSIS ECONOMICAMENTE 
ACONSEJABLE DE UN NUTRIENTE” 


Por MARIA ESTHER SUAREZ 


RESUMEN 


En este trabajo la autora, utilizando herramientas estadísticas y principios de la ciencia económica da 
una solución de lo que ““bautiza” como la “dosis económicamente aconsejable de un fertilizante o 
nutriente”. 

Destaca que la acotación que “efectúa al finalizar su investigación debe ser tenida profundamente en 
cuenta. 


SUMMARY 


In this work, with scientific rigor, but with a pragmatic objectiv, the author estimates has arrived to 
the solution of “the economical advisable dosis of a fertilizer”. She combines, to get her point, principles 
of these too important sciences: Statistics and Economics. 

The “marginal note” made at the end of this work, is importante taking into account. 


EXPERIMENTACION CON NUTRIENTES 
DETERMINACION DE LA DOSIS ECONOMICAMENTE ACONSEJABLE DE UN NUTRIENTE 


Este trabajo ha sido preparado cmo parte inicial de tan importante tema y al mismo 
tiempo de gran interés por tratarse de un problema factible de ser estudiado desde varios 
puntos de vista: económico, físico, matemático, estadístico y por supuesto agronómico. En 
este artículo nos introduciremos en el tema mediante un enfoque más quizás económico, 
técnico y estadístico matemático pensando en el gran interés que puede representar para el 
productor considerándolo desde el punto de vista de empresario y el problema que puede 
implicar para él la determinación de lo que podríamos llamar “la dosis económicamente 
aconsejable de un nutriente”. 

Nos hemos abocado a este tema especifico perteneciente al campo de la microeconomía 
por la relación estrecha que guarda, en lo que se refiere a implicancias económicas con la 
ciencia estadística. 


SUMARIO 
1. - Introducción. 2.— Las ecuaciones propuestas. El uso de la ley de Mitscherich en el 


análisis de experimentos con nutrientes. 3. — Fórmula para el cálculo de la dosis más 
económica de un nutriente. Conclusiones. 


lEscuela de Estadística, Facultad de Ciencias Económicas (Universidad Nacional de Rosario). Año 1974. 
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INTRODUCCION 


Todo análisis de experimentos con nutrientes, en particular cuando su estudio se hace 
desde un punto de vista económico y técnico, debe basarse en la hipótesis de que hay o 
existe una relación funcional entre las cantidades de nutrientes empleadas y el rendimiento 
obtenido de la variedad que-se trate. Pero qué es lo que sucede? , esta relación” no es 
conocida a priori pero puede ser aproximada en forma empírica con los datos numéricos que 
expresan el comportamiento de la producción relativa a las dosis de nutriente empleadas en 
la experiencia. 

Si expresamos con y, la producción agrícola, surge casi inmediatamente que debemos 
pensarla como una función de varios factores o variables, o sea que podemos darle una 
expresión simbólica escribiendo y = f (x¡, X2...Xn) en donde será la producción por unidad 
de área (medida por ejemplo en quintales de trigo o maíz por hectárea) y X;, X»,..X,, son las 
cantidades de los factores productivos empleados tales como abono, mano de obra, agua de 
riego, etc. 

A nosotros nos interesa encarar el estudio del comportamiento de la producción cuando 
hay un solo factor variable: el nutriente. Es decir, nos interesará determinar la ley analítica 
que liga la producción con las dosis diferentes de nutriente, por lo tanto lo expresaremos asi, 
buscamos y = f, (x,). 

Todo esto es del máximo interés para el experimentador y el productor cuando se plantean y 
se llevan a cabo experimentos con dosis crecientes de un nutriente. 

Antes de seguir con nuestro estudio, deseamos decir que si no se consiguen mantener 
constantes los demás factores en juego, el análisis estadístico revelará la importancia de ellos 
y los controlará a través de la técnica del Análisis de la Variancia. El investigador para 
determinar esta ley analítica y = f, (x,) toma como hipótesis de trabajo la llamada “ley de las 
proporciones variables”, llamada también “ley de los rendimientos decrecientes” o “ley de 
los crecimientos no porporcionales”, que en forma sencilla puede ser enunciada: “si la 
cantidad de un factor productivo fuera aumentada en incrementos iguales permaneciendo 
fijas las cantidades de otros factores, los crecimientos o incrementos de la producción 
correspondiente podrían crecer al comienzo mas decrecerían después de un cierto punto” 
(Spillman W.J. - 1924 - “The law of diminishing returns” - World Book Co. New York. — 
Heady E.O. - 1952 - “Economics of Agricultural Production and Resource Uses”. New 
Y ork.). 

Vamos a ilustrar esta ley por medio de la curva cuyas ordenadas representan el producto 
físico total y las abscisas el factor variable. Simbolizaremos a esta curva con PT (x) 
(deseamos aclarar que esta curva del producto físico total se la obtiene como P¿(x) dx, 
siendo P¿(x) la curva del producto marginal): 


y : PRA) 
la. Etapa | 22. Etapa | 3a. tapa 
i dy d 
Am : A iD E 
a él Po 


ao a SS —” 
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La curva PT (x) presenta una etapa inicial de rendimientos crecientes, esto es, PT (x) 
aumenta más que proporcionalmente ante la aplicación del factor. 

Viene luego una etapa de crecimiento de rendimientos decrecientes pero aún positivos, 
finalmente una última etapa de crecimientos de rendimientos nulos. 

Pero en realidad interesa conocer más el comportamiento de la curva PT(x), el del 
producto físico medio o sea Po (x) = e Pz (z) dz/z y el del producto físico marginal O sea, 
Pa (Xx), ya que del estudio del comportamiento de esta última surgirá la ecuación matemática 
desye 

El Po (x) Ó productividad media se refiere a la relación y/x, o sea, el producto físico 
total dividido por la cantidad del factor variable. Para todo punto de la curva PT (x), la 
productividad media está dada por la tangente del ángulo formado por el eje de abscisas y la 
recta que pasa por el origen. El producto físico marginal está dado por la derivada dy/dx. En 
la primera etapa, la curva P¿(x) es creciente, en la segunda etapa es decreciente. 

En los experimentos con nutrientes, casi nunca existe la primera etapa pues desde el 
comienzo los crecimientos de producción son en la gran mayoría de los casos, menos que 
proporcionales a las cantidades de nutrientes empleadas. Esto se debe sin duda al hecho de 
que las dosis del nutriente no son nunca muy pequeñas y ya el suelo posee una cantidad del 
abono en cuestión, lo que aparecerá luego como un parámetro en la ley analítica buscada. 
Este hecho apresura también el proceso de saturación de la producción. 


DETERMINACION DE LA RELACION FUNCIONAL ENTRE LA PRODUCCION 
Y LAS DOSIS DEL NUTRIENTE 


En el estudi o de los ensayos con nutrientes y con ayuda de los principios económicos 
que ya mencionamos, surge el problema de determinar a partir de los datos experimentales, 
la ecuación matemática de curva y = f, (x). Esto ha dado lugar a la proposición de 
diferentes ecuaciones y como generalmente sucede, ha sido elegida como “mejor” una de 
ellas. Sin embargo consideramos interesante mencionar algunas de las otras ecuaciones 
propuestas. 

En el problema de determinar la ley analítica nutriente-producción, es posible postular 

Varias ecuaciones como funciones de regresión de y en x y luego, llevado a cabo el 
correspondiente ensayo de hipótesis estadístico (con la correspondiente estadística P, cocien- 
te de dos variables), es factible encontrarse con el hecho de que ninguna hipótesis sea 
rechazada. Sin embargo, desde el punto de vista matemática; las ecuaciones pueden diferir 
enormemente como ya lo veremos. Esto explica uno de los aspectos característicos de las 
técnicas estadísticas. Un test de hipótesis puede conducir al rechazo de una hipótesis, pero 
nunca prueba o llega a probar que esta hipótesis sea la correcta, puesto que ya mencionamos 
que es siempre posible postular diferentes hipótesis que no son rechazadas por los test 
estadísticos. 
Lo crucial es que la elección entre las siguientes posibles hipótesis debe ser hecha en función 
del conocimiento específico del problema que da lugar a las observaciones. En el caso que 
nos ocupa, el razonamiento e hipótesis a utilizar para la determinación de y = f,(x) (que 
llamaremos simplemente a partir de ahora y = f (x)) será la siguiente: 

a) para una dosis de nutriente en una. parcela dada, la producción estará determinada 
principalmente por los siguientes factores: (a) características propias del nutriente; (6) la 
fertilidad natural del suelo y su contenido intrínseco de ese nutriente en cuestión y (y) la 
producción máxima teórica posible a obtener con valores crecientes de x. Veremos que sobre 
la consideración de estas hipótesis surge la ley analítica que dará la regresión teórica de y en 
x conocida como la ley de Mitscherlich, ley que cae dentro de las conocidas como “curvas 
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de crecimiento” y que caracterizan lo. fenómenos de crecimientos en el campo de la 
economía, física, química, población, etc. 

Es importante destacar nuevamente que la ley que va a representar la y = f(x) y que da 
la relación. funcional entre las dos variables x é y, se obtiene sobre la base de un 
conocimiento profesional del problema en discusión y esto es de fundamental importancia en 
lo que respecta a la fidedignidad de las extrapolaciones que se deseen hacer. 

Desarrollaremos entonces, primero en forma teórica en general, la deducción de una de 
estas curvas de crecimiento y veremos luego como este razonamiento lógico-deductivo puede 
aplicarse al caso concreto que nos incumbe. : 

Supongamos que la variable x simboliza el tiempo o la magnitud del factor que incide en 
el crecimiento de y en el fenómeno observado. Entonces, el coeficiente diferencial dy/dx 
simboliza la “tasa de crecimiento”, es decir, el crecimiento por unidad del tiempo (en 
nuestro caso por unidad de dosis de nutriente). 

Supondremos que el proceso de crecimiento puede ser caracterizado por la ecuación 
diferencial: 


dy/dx = f (x, y) 


la que indica que la tasa de crecimiento depende de ambos, de x, y del tamaño o cantidad 
obtenida de y. A partir de ahora trataremos sólo con el caso especial del tipo: 


dy/dx = h (y). g (x) (2.1) 
o bien: dy/h (y) = g (x) dx 
Integrando resulta: F (y) =G (x) 


y por medio de esta ecuación, y, queda determinada como una función de x. 
Introducimos ahora en (2.1) el supuesto que: 


dy/dx =(A - y) g (x),(0< y<M) (2.2) 


Por qué hemos introducido esta h (y) expresada en esta forma? . Ya lo dijimos antes, esta 
ecuación que incluye el parámetro A como cota superior de y, indica que el crecimiento está 
limitado ya que A representa el valor externo que podría tomar y. 

Ella indica que para un tiempo dado, la tasa de crecimiento es proporcional al “tamaño 
remanente”, esto es, al tamaño máximo menos el tamaño alcanzado y una función del 
tiempo. Si ahora consideramos el “coeficiente diferencial logarítmico”: 


dilogiyaiza le dy 
dx y dx 


que representa “la tasa de crecimiento relativa”, es decir, el incremento proporcional por 
unidad del tiempo, nuestra ecuación vendrá: 


d ae A O) les) 


dx (MY) 


_dlog(A— y)= g (x) 
dx 
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y la solución general de nuestra ecuación diferencial será: 
IN la do a 6 (a) A 
A e Ela , y O A, 


pero ahora, por las condiciones del problema, la constante indeterminada del problema toma 
en este caso particular el valor H¡ = A, pues cuando x = 0, y = 0 y esto implica que H, =A, 
de donde resulta que: 


A Ie On (2.3) 


Falta ahora explicitar a g(x). Introduciendo funciones especiales para g(x), por ejemplo, 
podrían ser: 


En A AA Ó g2 (Xx )=Bo + Brlx + Ba/x” 


y obtendríamos así un número de ejemplos de curvas de crecimiento. Para nuestro caso 
vamos a admitir el más simple, que g(x) = f, obteniendo entonces como solución: 


Alo e] MES (2.4) 


Para B > 0 ésta es una función creciente de x, teniendo para x>” a A como asintota. 

Deseamos destacar que desde el punto de vista estadístico surge un problema especial en 
lo que se refiere a la estimación de los parámetros de esta ecuación (A y y f) ya que las 
observaciones generalmente no son estadísticamente independientes, entonces los supuestos 
que dan las bases del Analisis de la regresión generalmente no se mantienen. Los métodos 
estadísticos a ser-usados para el tratamiento de datos provenientes de procesos de crecimien- 
tos, aún no están bien clasificados. 

Elementos sumamente valiosos pueden obtenerse en la bibliografía dada al final. 

En estos ensayos con nutrientes, que presentan todas las características de las curvas de 
crecimientos, hay entonces motivos poderosos para evitar el uso de regresiones parabólicas, 
dándose preferencia a una regresión de y, en x, que con razonamiento análogo al general 
presentado en páginas anteriores, conduce a la llamada ley de Mitscherlich, cuya expresión 
analítica es: 


AO 


en la cual: A es la producción máxima teórica posible para valores muy elevados de x, (dosis 
de nutriente), c, llamado coeficiente de eficacia, esun parámetro típico de nutriente en 
cuestión y b, corresponde a la fertilidad natural del suelg, sin el aporte del abono en estudio. 
En los países en los que rige el sistema métrico, en la ecuación de Mitscherlich x, se 
expresa en quintales métricos por ha.; c, es entonces expresado en ha./quintal métrico. 
Para hacer un poco de historia (lo que siempre resulta ameno o instructivo) podemos 
decir, que los ingleses prefieren otro tipo como'la ecuación nutriente-producción. 
Algunos como Stevens, prefieren dar a la ley de Mitscherlich la forma: y =a + f£* o lo que 
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es lo mismo, y = m - ar*, ésta debida a Spillman y que corresponde exactamente a la 
ecuación propuesta por Mitscherlich. 

Como vemos, la ecuación de Mitscherlich admite una asíntota horizontal, de ordenada 
igual a A. Esto está perfectamente de acuerdo con lo que sucede en la práctica en el uso con 
nutrientes, pues éstos a partir de una cierta dosis en adelante llegan a producir un efecto de 
saturación. (Esto es de sumo interés para entrar luego a considerar el punto de la determina- 
ción de la dosis económicamente aconsejable de un nutriente). 

Aunque mencionamos quepara estos casos, el uso de polinomios no conviene para el 
análisis de estos experimentos (considerando los hechos reales que nos llevan a la formula- 
ción de hipótesis concordantes con la realidad), diremos que otros autores, como Heady por 
ejemplo, en “Economic Production Functions”, prefieren utilizar un trinomio de segundo 
grado: a + bx +b2x?. ! 

Esta ecuación conduce a cálculos generalmente más simples pero tiene defectos graves 
que hacen desaconsejable su uso, como ser: 
l)es completamente impropia para una extrapolación, lo qe no sucede con la ecuación de 
Mitscherlich. 2) Los parámetros con excepción de a, no tienen una interpretación objetiva 
conveniente. 

Al revés de esto, la ecuación de Mitscherlich aunque presenta las dificultades en las 
estimaciones de sus parámetros, éstas son compensadas por ventajas que aparecen posterior- 
mente y que la recomiendan para su uso en forma ya más general. Ella exige, para su 
aplicación conveniente, experimentos de gran precisión o mejor aún grupos de experimentos 
numerosos. En estos casos, los valores de c, obtenidos para cada nutriente son practicamente 
constantes. Nota: Es desaconsejable su uso, cuando el suelo reacciona muy mal ante el 
nutriente en estudio. 

Según mencionan en sus artículos los Profesores A.C. Zaggatto y F. Pimentel Gómes, en 
Experimentos de campo en regiones tropicales y subtropicales permiten tomar como muy 
razonables las estimaciones de c, siguientes: c = 0,49 ha./quintal métrico para el nitrógeno; c 
= 0,88 ha./quintal métrico para el fósforo y la potasa (K,0) y para el estiércol c = 0,018 
ha./ton. 

(Hasta el momento de realizada esta investigación no se ha posido constatar estadísticamente 
el valor estimado de los parámetros por falta de información solicitada y tiempo. Será una 
deuda para el inmediato futuro). 

Para mayor facilidad en la aplicación de la ley de Mitscherlich (lo mismo que en el caso 
de los polinomios ortogonales), se deben adoptar niveles de nutriente igualmente espaciados. 

La gráfica de la ecuación de Mitscherlich es: 
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LA FORMULA PARA EL CALCULO DE LA 
DOSIS ECONOMICAMENTE ACONSEJABLE 
DE UN NUTRIENTE 


Este problema de la determinación dela dosis “Óptima” de un nutriente es sumamente 
importante y surge como una consecuencia inmediata e la ley deMitscherlich, representativa 
de la regresión producto-nutriente. Mencionamos ya anteriormente que esta curva refleja la 
existencia de un proceso de saturación de la producción ante la aplicación de dosis crecientes 
del nutriente, por lo tanto es casi evidente el deseo de determinar un valor “óptimo” de la 
variable independiente x (o sea la dosis de nutriente) para evitar entrar en esta etapa de 
saturación, esto es importante aún prescindiendo del problema monetario, el que sabemos 
tiene una gran importancia real, por lo tanto razonaremos este problema basándonos en los 
conceptos previos pero incluyendo también el problema económico, para poder llegar así al 
valor de x, que será el valor buscado. Veremos que mediante aproximaciones, este valor 
resultará igual a: 


E E a log. A (3120) 
De 6 t = Xu 


donde u es el aumento de producción del producto agrícola considerado obtenido con la 
dosis X,, de nutriente; v, es el precio unitario de este producto pagado por el productor, t, es 
el precio unitario del nutriente y c, es el coeficiente de eficacia. 

Es evidente que siendo P, el precio unitario del producto y P, el precio unitario del 
nutriente, sólo continuaremos aplicando dosis adicionales de nutriente cuando el crecimiento 
o incremento Py dy en el valor de la producción, no son inferiores al incremento P_dx, es 
decir, siempre y cuando el rendimiento marginal sea mayor que el costo marginal: 


Pod Pd dy dx =P Pp, (352) 


En los estudios con nutrientes dy/dx es una función decreciente de x, de donde surge 
que siendo x* la solución de la ecuación: 


dy/dx =P,/Py (37. Si) 


debemos usar una dosis x < x*, esto es, x* es la dosis máxima económicamente conveniente, 
para cualquier dosis de nutriente inferior a x* será también lucrativa. 

Antes de resolver esta ecuación (3.3) vamos a considerar algo que llamaremos “un factor 
de confianza o corrección”, sea f>1 y teniendo en cuenta los efectos residuales de 
nutrientes en años sucesivos, obtenemos la ecuación: 


Pp, e a de de Pd, 


donde y; (i=1, 2, 3, 4, ...) se refiere a la producción en el año iésimo o a la i-ésima 
cosecha. Supuesta adoptada la ley de Mitscherlich, encontramos mediante algunas hipótesis 
simplifacoras del problema, la siguiente ecuación: 
1 l— : 
O A o pp A e] (3.4) 
c f (f—h) P, log e 
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donde fh mide el efecto residual del nutriente. En el caso que f — h= 1, podemos deducir la 
fórmula aproximada: 


A P, C 


q O (55) 
E LL log e 


J 


a [ 
E 


la que para f= 1 corresponde exactamente a la solución de la ecuación (3.3 ). 

En el caso general de tener dos niveles de fertilizante (cero y x,,), con incremento de 
producción, u, causado por el nutriente, esta fórmula (3.5) puede ser aproximada bastante 
bien por: 

E E la (3.6) 
2 C EX 


AS 


Donde (3.6) es igual a la primera fórmula dada en (3.1). 

Pero para que se pueda aplicar la dosis “Óptima” y máxima x*, será necesario que no 
haya limitaciones de capital a disposición del productor, esto es, que el capital sea suficiente 
para la compra de esa cantidad óptima. Supongamos por ejemplo, que el productor sólo 
dispone de un monto o capital C a ser empleado en la compra del nutriente. Es decir, que 
tenemos que introducir una restricción adicional a la (3.2), que es: 


Siendo x** la dosis de nutriente a ser empleada por el productor, estará limitada por las 
desigualdades x < x*, x < x**, luego no podrá exceder a los dos valores x* y x**, 
CONCLUSION 
Estas leyes que me dan la producción según el nutriente utilizado, tienen el valor de las 
leyes económicas cuando se introduce la referencia de ceteris paribus, lo que quiere decir que 


las leyes que se dan son paramétricas, es decir, están condicionadas a la naturaleza y medida 
de las otras variables, observación ésta que no siempre se tiene en cuenta. 
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APLICACION DE LA TEORIA DE VARIABLE COMPLEJA AL CALCULO DE 
AUTOVALORES EN GUIAS DE ONDA DE SECCION ARBITRARIA 


Por PATRICIO A. A. LAURA* 


RESUMEN 


El uso de guías de onda de sección no convencional se ha hecho común en la última década. Diversos 
métodos analíticos y numéricos han sido utilizados para determinar modos de propagación, autovalores, 
etc. Un método que ha sido usado por varios investigadores en Europa y Estados Unidos de Norte 
América es el de transformación conforme, por el cual se convierte a la sección en una más simple y luego 
se aplica alguna de las técnicas de residuos ponderados para resolver el sistema diferencial transformado. 
En este trabajo se presentan algunos resultados de investigaciones recientes realizadas por el autor y sus 
colaboradores para secciones simples y doblemente conexas. 


INTRODUCCION 


Las guías de onda de sección circular y rectangular no siempre ofrecen una solución 
óptima desde el punto de vista practico. Las investigaciones tanto experimentales como 
analíticas sobre secciones de forma compleja son numerosas. 

Los métodos analíticos que han sido utilizados por diversos investigadores son*: 


(a) transformación conforme 

(b) diferencias finitas 

(c) colocación por puntos en el contomo 
(d) método de Rayleigh-Ritz 

(e) método de Galerkin 

(f) método de momentos 

(g) técnicas de perturbación 

(h) métodos de impedancia 

(1) elementos finitos. 

En el presente trabajo se analiza el cálculo de autovalores en guías de onda de sección com- 
pleja, mediante el método de transformación conforme y técnicas variacionales. Este enfoque 
ha sido utilizado por el autor y sus colaboradores en la solución de varios problemas de deter- 
minación de autovalores en guías de onda de sección “no circular” y “no rectangular”. 

Es importante mencionar los trabajos de Meinke2,3 y Tischer%-7, quienes han hecho uso 
del método de transformación conforme en el estudio de propagación de ondas electromag- 
néticas en varios tipos de guías de onda. 

Ambos autores3P transforman a la sección transversal de la guía de onda en un 
rectángulo siendo posible calcular frecuencias críticas, campos, corrientes; potencia transpor- 
tada y atenuación de los diversos modos de propagación. El método ha sido desarrollado en 
forma independiente por ambos investigadores consistiendo la diferencia fundamental entre 
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los enfoques en las técnicas de solución de la ecuación transformada de Helmholtz. 

Meinke y Tischer mencionan la aplicabilidad de procedimientos que hacen uso de la 
técnica del tanque electrolítico para obtener la transformación deseada. 

En Ref. 7, Tischer utiliza la transformación de Schwartz - Christoffel para resolver un 
importante problema de la teoría de microondas. Bava y Perona* resuelven un interesante 
problema de microondas siguiendo el método de Tischer”, pero utilizando tres transforma- 
ciones sucesivas. 

El autor y sus colaboradores?-13 han usado el método de representación conforme, 
transformando a la sección dada en un círculo unitario en el caso de regiones simplemente 
conexas y en un anillo circular si se trata de regiones doblemente conexas. Resulta así un 
sistema diferencial en el cual las condiciones de borde pueden ser satisfechas idénticamente, 
pero donde la ecuación diferencial debe ser resuelta en forma aproximada. Esto es llevado a 
cabo mediante el método de residuos ponderados. 


TEORIA 


En el caso de una guía de onda prismática perfectamente conductora, las ecuaciones de 
Maxwell se reducen a la ecuación de Helmholtz: 


YY +k?%w =0 (1) 
donde “es una función escalar k= KG Koko. Es > o 20 o ORIO 
de onda: V? es el operador Laplaciano en dos dimensiones y w es la frecuencia circular del 
modo de propagación. 

Las condiciones de borde que debe satisfacer Y son: 


y [L (xy) = 0] =0 (modos transverso magnéticos) 24) 
NY LO 

Ay) 00 (modos transverso eléctricos) (2b) 
ón 


donde L (x,y) = O es la relación funcional que define el contorno de la sección en el plano 
w (w=x + iy) y n es la normal al borde. 
Es conveniente escribir a la ecuación (1) en la forma: 


2alr . 
A (3) 
gw dWw 
donde w = x - iy. 
Sea 
w=f(E) (4) 


la función analítica que transforma a la región (simple o doblemente conexa) es una más 
simple en el plano ¿ (círculo unitario en el caso de una región simplemente conexa y anillo 
circular si se trata de una región doblemente conexa). De las ecuaciones (3) y (4) se obtiene: 
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2 
a e 1 (E)1? =0 (s) 
08 0É 
Las condiciones de contorno transformadas son, en el caso de una región doblemente 
conexa: 


Y(5E Jl =00 == 19) (6a) 


A 0% 
dr Mene) ly 0 (LD) (6b) 


En este trabajo solo se consideran modos transversos magnéticos de propagación. El caso 
de modos transversos eléctricos, puede ser analizado también siguiendo el procedimiento de 
HinaY. 


SOLUCION DEL SISTEMA DIFERENCIAL TRANSFORMADO 


a. Metodo de Galerkin 


En el método de Galerkin se aproxima la solución de (5) mediante una sumatoria: 


N 
2 On Un (EE) (7) 
ne 
donde las Y, forman un conjunto completo de autofunciones que satisfacen idénticamente 
las condiciones de contorno, Substituyendo (7) en (5) resulta una expresión residual o función 
error € (£, E), que en general será diferente de cero. El método de Galerkin requiere que e (£, 
£) sea ortogonal con respecto a cada una de las funciones coordenadas Wy (£, E) en el, 
dominio transformado D: 


E (E EDI 0 (8) 


Resulta así un sistema homogéneo de N ecuaciones lineales en los coeficientes b,,. Para que 
este sistema tenga una solución no trivial, es necesario y suficiente que el determinante de 
los coeficientes de las incógnitas sea idénticamente nulo, obteniéndose así una ecuación-de- 
terminante en los autovalores k. 

Sea la configuración que se muestra en la Figura 1, cuya transformación en el plano W 
está dada por la función analítica; 


AE am o E == 1010 (9) 


donde m = 1/6. 
Tomando a (7) en la forma: 
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Y=E bp11-(5.Fp (10) 


y reemplazando (9) y (10) en la ecuación diferencial gobernante (5) se obtiene la correspon- 
diente expresión de la función error e(r,0). Tomando n= N = 1 resulta: 


do 
a 


Ko 1 


y haciendo N = 2 resulta?? : 
2.388 
Kko1 IS 

a 


Cabe mencionar el hecho de que este valor debe ser considerado como una primera 
aproximación, ya que la ecuación (10) contiene únicamente a la variable radial. 


b. Método de Colocación a lo Largo de Arcos 


En este método se requiere que la integral de la función error a lo largo de N arcos 
concéntricos se anule: 


ade E) 2. 1000 0= 12) (11) 


resultando así un determinante-ecuación en los autovalores k. 
En Ref. 16 se presentan resultados tomando: 


N a 
US 2, Bao Lan (E.E)?] (12) 


donde J, denota la función de Bessel de primera clase y orden cero y los parámetros Q, son 
sus raíces. 

Tomando "N= Py er, 0.0230. 305075 resulta! : 
ARO, 

a 


Ko 
c. Método de Colocación en Subdominios 


El criterio de minimización del error consiste ahora en requerir que la integral de la 
función error en N subdominios se anule: : 


E E) 0 = 0 (e B22N) (13) 


Tomando N = 2, el valor obtenido en la Ref. 16 es k,¡ = 2.394. 
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d. Determinación de Límites Superior e Inferior por el Método de Kohn-Kato 


Este método es de uso común en física nuclear y fue extendido a problemas bidimen- 
sionales en Ref. 17. La aplicación del método es sumamente laboriosa obteniéndose en defini- 
tiva los siguientes valores!” : 


2.384 
ad 


2-399 


< Ko < a 


Desde el punto de vista práctico, los límites calculados son tan cercanus que el corocir. iento 
de los mismos equivale a conocer el valor exacto. Es importante mencionar, el hecho de que es:-. 
es la única configuración compleja para la cual se conoce su primer autovalor con tanta pre- 
'cisión. 

Comparando los resultados obtenidos (Tabla 1), se concluye que todos los valores con- 
cuerdan con muy buena precisión. El método de Galerkin da un valor muy cercano a los lími- 
tes obtenidos mediante la técnica de Kohn-Kato. 


OTRAS APLICACIONES 
Sea ahora el caso de una guía de onda de sección cuadrada para la cual se conocen los 


autovalores exactos. 
La función transformación está dada por la integral elíptica: 


É 
w=-1082 So (1+ 42 dt (14) 


donde a es la apotema del cuadrado. 
Expandiendo a (14) en series de Mac Laurin se obtiene: 


Z 
! 


a(1.08£—0.108£% + 0.0455? — 0.026£* +...) (15) 


Tomando N = 2 en la ecuación (10) se obtiene mediante el método de Galerkin: 


ZA 5.176 
a a 
a a 


y tomando N = 5 en (12) se obtiene por colocación a lo largo de arcos: 


23221 x 5.069 
Ko = AN , ko = a 


siendo los valores exactos: 


2.221 4.967 
Ko = a , Ko = = 
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La Tabla 2 muestra resultados obtenidos para diferentes secciones poligonales regulares'' 
Los resultados obtenidos para el primer autovalor son prácticamente los mismos. En el caso 
del segundo autovalor ko , el valor obtenido mediante la técnica de colocación a lo largo de 
arcos es 2% más alto que el valor exacto para la sección cuadrada, mientras que el método 
de Galerkin da un valor 4%-más alto. 

La Tabla 3 muestra una comparación de resultados obtenidos por diferentes autores para 
una configuración exagonal. Se observa que hay una variación drástica de los autovalores con 
la densidad de la malla que se utilice en el método de diferencias finitas. 

En el caso de una región doblemente conexa (Fig. 2), Referencias 12 y 13, describen 
investigaciones analítico-numéricas sobre regiones poligonales con orificios circulares. 

Los autovalores Ko, y Kg han sido calculados utilizando la siguiente expresión polinó- 
mica: 


Y (1,0) = Ar AU=D(E-A) + Az (1-12)(-A) (16) 


donde el radio externo rz es tomado igual a la unidad y el radio interno r, es/igual a A. 
La Figura 3 muestra la variación de los autovalores Ko, y Kg, en función del parámetro 
apotema/radio del contorno interno para la configuración de la Figura 2. 


CONCLUSIONES 


El método de variable compleja da valores numéricos muy precisos y ofrece excelentes 
condiciones de estabilidad. El cálculo de modos superiores de propagación es también 
posible, pero se torna más complejo. | 

Determinación de la función analítica que efectúa la transformación, es en general, una 
tarea difícil pero cabe mencionarse el hecho de que existan numerosos métodos numéricos 
para hallarla. 
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TABLA | 
Comparación de resultados obtenidos mediante diferentes métodos (Figura 1) 


Método utilizado . Ko¡.a Autor 

Método de Galerkin 2.388 Laura y Shahady?* 
Técnica de Colocación Baltrukonis, Chi, Casey 
a lo Largo de Arcos DS y Laura!* 

Técnica de Colocación Baltrukonis, Chi, Casey 
en Subdominios 2.394 y Laura!* 

Método de Kohn-Kato 2.384< x < 2.385 Baltrukoni, Chi y 


Laura!” 


Sección 


Cuadrada 


Pentagonal 


Exagonal 


Eptagonal 


Octogonal 


Circular 


MODO 


Kg .a 


APLICACION DE LA TEORIA DE VARIABLE COMPLEJA. .. 


TABLA Il 
Autovalores calculados para diferentes secciones poligonales regulares. (Laura y Faulstich, Ref. 11) 
Autovalores 
loci - Colocación a lo Método de 

Largo de Arcos Galerkin 
ko, .2 ZA Z 12-221 ZED ZA 
ko. .a 4.967 5.069 5.176 
ko, .a a 2.285 2.285 
koz .a AENA 3.239 5.19 
ko, .a O EAS 1 NA 2.318 
ko, .a o 5.313 5.278 
ko¡ .a EA 2.339 2.339 
ko, .a A 5.367 EE 
ko1 .a ES 4 UR 2.339 2.353 
Ko2 .a E 5.405 ————— 
ko: .a 2.4048 2.4048 2.4049 
ko. .a 5.5201 5.5201 3.321 

TABLA 1!I 
Comparación de autovalores de una guía de onda de sección exagonal 
Colocación 
Diferencias a por Puntos Colocación a E Galerkini1 
Número de Celdas en el Largo de Arcos 
Contorno12 
1) 19 6 12524 
DTD DDD 2.313 2.317 2.318 
3.90 4.62 — 5.54 4.81 — 5.313 5.278 


Ko .a 


Koa .a 


5.85 8.20 — 8.33 = 
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FIGURA 2 - TRANSFORMACION APROXIMADA DE 
UNA REGION DOBLEMENTE CONEXA 
(CUADRADO CON ORIFICIO CIRCULAR) 


(Refs. 12 y 13). 


APLICACION DE LA TEORIA DE VARIABLE COMPLEJA. . . 
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TRICOS EN EL PLANO 8) DE LA 
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EL RIESGO SISMICO EN EL TERRITORIO ARGENTINO 


Por FERNANDO VOLPONI!' 


RESUMEN 


Es frecuente considerar que en las zonas que se sabe que son sísmicas deben ser adoptadas medidas 
de prevención, pero que en las conocidas o supuestas como nó sísmicas no es necesario hacer lo mismo. 
La clasificación de las zonas en sísmicas y no sísmicas se basa generalmente en observaciones de carácter 
estadístico. Pero la información disponible actualmente no es abundante especialmente para las regiones 
poco activas, y se suele considerar como no sísmicas zonas que en realidad sí lo son en cierto grado. 
Vemos que de cuando en cuando se producen terremotos en lugares donde no se esperaba que pudiesen 
ocurrir. 

Esto nos hace pensar que es importante investigar la sismicidad no solamente en las zonas activas si 
no también, y con la misma intensidad, en las zonas poco activas del territorio nacional. 


SUMMAR Y 


It is well known we have to take some preventions against earthquakes in sismic zones, and, 
frequently, it is supposed that it is not necessary to do the same in non seismic ones. 

But the classification of zones in seismic and non seismic ones is based on statistical observations, 
and we do not have enough data to do that correctly, especially in relation with the less active regions. 
Frequently active zones are considered as non active ones and from time to time, we have earthquakes in 
places were nobody expected they could occur. 

This means that it is necessary to investigate the less active zones with the same deep interest as the 
active ones. 


INTRODUCCION 


El problema sísmico argentino no está solamente en la Cordillera. Para enforcar el 
problema en su aspecto general hay que hacer algunas consideraciones previas que no son de 
carácter puramente geofísico. 

El seguro contra algún siniestro depende de dos factores. Uno es la probabilidad de 
ocurrencia del siniestro, y el otro es el valor que está expuesto al riesgo. 

En el caso de terremotos, el primer factor es un índice numérico mediante el cual 
representamos la sismicidad de una zona, y el segundo es un valor que depende del grado de 
desarrollo de la misma zona. Desde el punto de vista de los intereses en juego, puede suceder 
que resulte más importante investigar el problema en una zona poco activa pero densamente 
poblada que en otra más activa pero escasamente poblada. Estas consideraciones alcanzan 
una gran aplicación en el ámbito del territorio nacional, que se extiende de Este a Oeste 
comprendiendo grandes áreas, a lo largo de las cuales varía fundamentalmente tanto el índice 
de sismicidad, como la densidad de la población o su grado de desarrollo. 


l Investigador del Consejo Nacional de Investigaciones Científicas y Técnicas. Director del Instituto 
Sismológico “Zonda” de la Universidad Nacional de Cuyo. Profesor titular de GEOFISICA APLICADA . 
en la Facultad de Ingeniería de la misma Universidad. 
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ESTIMACION DEL INDICE DE SISMICIDAD 


Es difícil estimar correctamente la probabilidad de ocurrencia de un terremoto. Por 
ahora, no hay otra alternativa que conformarse con un resultado al cual hay que asignar un 
grado bastante elevado de incertidumbre. 

Tenemos la posibilidad de recurrir a tres criterios diferentes, y finalmente hacer una 
integración de los resultados finales de cada uno de ellos. 

a) El primero se basa en la interpretación estadística de todos los terremotos que se han 
registrado en la zona en estudio. 

b) El segundo se basa sobre la observación de determinadas estructuras geológicas o morfoló- 
gicas, que muestran la existencia de condiciones de inestabilidad tectónica. 

c) El tercero se basa sobre datos geofísicos que también demuestran la presencia de condicio- 
nes de inestabilidad. 


a) Mediante el análisis estadístico de todos los terremotos registrados en una zona, se 
puede calcular un índice de sismicidad correspondiente a la misma. Cuanto más abundante es 
la información, tanto más los resultados se aproximarán a los valores deseados. 

Pero la Sismología es una ciencia moderna, y los datos disponibles suelen ser escasos. Se 
distribuyen en un intervalo de tiempo demasiado corto como para definir, con precisión, la 
curva de distribución de los acontecimientos. 

Por otra parte, no existe uniformidad zonal en la recolección de los datos. Algunas áreas 
tienen una buena densidad de sismógrafos, pero otras se encuentran demasiado alejadas de 
cualquier equipo sismográfico. 

Los mapas de sismicidad global que ahora conocemos, corresponden a una sismicidad de 
terremotos fuertes (magnitud mayor que 5 Richter). Mapas globales correspondientes a la 
actividad de terremotos débiles todavía no se conocen. Estos ampliarian notablemente el 
conocimiento de la actual Geodinámica del planeta. 


b)la existencia de fallas activas, y la comprobación de movimientos recientes en las 
capas superficiales de una zona constituyen evidencias geológicas de la presencia de condicio- 
nes de inestabilidad en la misma. | 

La evidencia geológica más importante, que demuestra la presencia de condiciones de 
inestabilidad, y por consiguiente de actividad sísmica es, indudablemente, el borde continen- 
tal representado morfológicamente, en América del Sud, por la Cordillera de los Andes. Toda 
la Cordillera es sísmicamente muy activa, pero la sismicidad vá disminuyendo hacia el 
interior del continente. 

Si la condición de borde es la causa que determina la sismicidad de la Cordillera, se 
puede pensar que, si actualmente existe o se observa menor actividad en una parte de la 
Cordillera que en otras, tal hecho debe considerarse como una situación transitoria. Y que 
una zona hoy tranquila, puede ser mañana muy activa. 


c) Entre: los datos. geofísicos útiles en este problema están los que nos puede suministrar 
la Gravimetría. 


La presencia de anomalías isostáticas indica la existencia de una irregular distribución de 
masas, O sea, la falta de equilibrio isostático. La falta de equilibrio isostático presupone la 
- existencia de estados de tensión en las capas sólidas de la litósfera. Y dicho estado de tensión 
puede dar lugar a movimientos relativos entre bloques y así originar ondas sísmicas. 
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LO QUE SE CONOCE AHORA EN EL TERRITORIO ARGENTINO 


a) Del mapa de epicentros 


Si observamos el mapa de la Sismicidad Mundial, p.e. el publicado por la “COAST AND 
GEODETIC SURVEY” (1961-1969)! , vemos que en el continente Sudamericano la actividad 
sísmica se produce a lo largo de una franja que sigue la costa del Océano Pácifico, que 
comprende justamente la Cordillera de Los Andes; su ancho es variable, y alcanza su valor 
máximo a la altura de la Puna de Atacama, donde llega a unos 800 kilómetros. 

Al Este de dicha franja la actividad va decreciendo paulatinamente. 

Si observamos, en el mismo mapa, la distribución de la actividad en la parte norteameri- 
cana del continente, especialmente en el territorio de los Estados Unidos, vemos un panora- 
ma diferente. Vemos que la actividad sísmica se distribuye mucho más uniformemente a lo 
largo de todo el territorio mostrando, además, algunas concentraciones de actividad tanto al 
Este como al Oeste. Sabemos que las condiciones geodinámicas existentes en la parte Norte 
son diferentes de las existentes en el Sur; pero quiero hacer notar que las diferencias que se 
observan se deben, en gran parte, al hecho de que, mientras en los Estados Unidos abundan 
las estaciones sismológicas tanto al Este como al Oeste, en Sudamérica no sucede lo mismo. 
Existe un buen número de estaciones al Oeste (en Chile, Perú, Bolivia, Ecuador, Colombia y 
Argentina), pero hay muy pocas sobre la costa atlántica. El número de estaciones que se 
encuentran instaladas sobre las zonas del Este de América del Sud, son insuficientes para 
obtener la información que hoy se usa para calcular los sismos que en ellas se producen. Por 
consiguiente dichos sismos no figuran en los catálogos universales y tampoco resultan 
representados en los mapas de la sismicidad global. Mientras que en los Estados Unidos un 
sismo de magnitud 3 Richter aparece representado en los mapas, en Sud América no 
aparecen los de magnitud hasta 6 Richter si se producen en las regiones orientales. 

Por lo que se acaba de decir, en la Argentina es muy poco lo que se conoce de las 
características sísmicas del territorio, si se exceptúa la zona cordillerana. Prueba de ello es 
justamente que cada vez que se produce un terremoto un poco fuerte, en algún lugar de las 
llanuras orientales, causa enorme sorpresa y estupor. 

En resumen, del mapa de la actividad sísmica del territorio argentino obtenemos ahora 
una información ya útil pero todavía demasiado incompleta. Se puede decir que ya es 
suficiente la que corresponde a las zonas muy activas, pero decididamente insuficiente la que 
corresponde a las áreas poco activas. 

En un intento para mejorar la información disponible, se han investigado, por vía 
bibliográfica, los terremotos antiguos que se produjeron en el territorio argentino y zonas 
vecinas. Los más importantes son los que se presentan sintéticamente en la TABLA Il. 
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TABLA | 
Lista de Terremotos ¡votables que se produjeron en el 
Territorio Argentino antes de 1900 
FECHA ZONA OBSERVACIONES 
31 Oct 1527 Cerca de San Nicolás Relato de la Expedición de Caboto (4). 
1692 Esteco (Salta) Desaparece parte de la población. (4). 
18 Oct 1844 Salta, Jujuy, Tucumán Víctimas y Daños (8). 
y Sgo. del Estero 
4 Oct 1845 Lobos, Navarro, Chivil- Informe del Dr. Francisco Javier Muñiz “Gaceta 
coy (Prov. Bs. As.) Mercantil” N“ 6716 del 26 Feb. 1846, Buenos 
Aires. 
9 Ago1848 Montevideo y zonas Fué seguido de réplicas importantes. (8). 
vecinas 
1849 San Luis (8) Terremoto fuerte (8). 
16 Nov 1857 Mendoza Terremoto fuerte con dos réplicas importantes 
(8). 
Jul 1858 Salta Terremoto fuerte (8). 
20 Mar 1861 Mendoza Destrucción de Mendoza (10.000 víctimas, 
muerte del naturalista Bravard) (8). 
11 Nov 1876 San Juan y Mendoza Derrumbe de edificios (8). 
27 Mar 1877 Córdoba Alarma en la ciudad (8). 
12 Feb 1880 San Luis Se hundió un cerro “*? ” (8). 
17 Ago 1880 Mendoza Muertos y Daños (8). 
21 Ago 1880 Mendoza Muertos y Daños (8). 
5 Jun 1888 La Plata hasta Campana Alarma en Nuenos Aires, Informe del Dr..F. 
(Prov. de Bs. As.) Ameghino, “La Nación” del 13 de Junio de 
1888, Buenos Aires. 
27 Oct 1894 San Juan, La Rioja Gran Terremoto Informe del Dr. Guillermo 


Bondembender presentado al Rector de la Univ. 
de Córdoba, Dr. Telasco Castellanos. 


b. De los Datos Geológicos 


El terremoto del 15 de Enero de 1944 que destruyó la ciudad de San Juan provocó un 
deslizamiento en un plano de falla que se pudo observar en la superficie. Medido, alcanzó un 
valor de unos 35 centímetros en la dirección de la vertical. Es una falla inversa que tiene 
dirección NE-SO, y se ha levantado el bloque SE. Desde entonces esa falla recibe el nombre 
de falla de Las Lajas. Posteriores investigaciones del lugar me han permitido descubrir que se 
trata de una falla importante preexistente. El hecho que en esa zona, en la superficie, existen 
materiales clásicos modermos, ha impedido que la falla se observara con anterioridad. El 
rechazo de la falla, que actualmente se puede medir con solo comparar los desniveles de los 
bloques a ambos lados de la misma, podría alcanzar unos 8 metros. La traza de la falla se 
hace muy evidente en vistas aéreas, pero también se puede identificar desde la superficie. 

Una falla activa, o una falla moderna, constituye una manifestación geológica de 
actividad sísmica. 

Se supone, con bastante fundamento, que las llanuras orientales del territorio nacional, 
se encuentran divididas en grandes bloques, separados por sistemas de fallas. Si bien muchas 
de estas fallas han sido estudiadas aprovechando sus afloramientos, es indudable que otras 
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permanecen cubiertas por sedimentos recientes, y de las cuales poco se puede saber. 

Los importantes terremotos de Sampacho del año 1934 y sus precursores, indican 
claramente la presencia de una gran falla activa de direción Norte-Sur aproximadamente. 

Tal vez suceda algo semejante con los terremotos de Santa Rosa (La Pampa), de Arias 
(Córdoba), de Tostado (Sta. Fé) y de Las Breñas (Chaco) véase la Tabla II, y figura 1. 

Mas al Este tenemos el importante terremoto del día 21 de Enero de 1948 cuyo foco 
debe haber estado entre las localidades de Goya y Esquinas sobre el curso del rio Paraná y 
cuya magnitud debe haber sido de alrededor de 6 Richter. Fue sentido en una zona muy 
extensa de la Argentina y del Uruguay, y causó alguna alarma en la ciudad de Buenos Aires. 

Si bien todavía es poco el peso que se le puede asignar a los factores geológicos en el 
problema de determinar el índice de sismicidad de las llanuras orientales del territorio 
nacional, ellos constituyen igualmente elementos de juicio muy importantes. 


c) De los datos geofísicos 


En correspondencia con las zonas más sísmicas de la tierra existen fuertes anomalías 
gravimétricas. Parecería que las fuerzas que tratan de restablecer el equilibrio isostático 
fueran las responsables (por lo menos en parte) de los movimientos sísmicos. 

Una de las partes en las cuales se observa más claramente este fenómeno es en la costa 
occidental de Sud América, o sea, en el borde continental representado por la Cordillera de 
Los Andes, y en especial a la altura de la Puna de Atacama. 

Por el análisis del mecanismo del foco de los terremotos, se infiere que sobre los bloques 
de la litosfera actúan fuertes presiones en dirección horizontal. Estas presiones hacen que 
resulte difícil la interpretación del desequilibrio isostático. 

Las características de tales presiones son casi desconocidas, no se sabe como varían 
según el azimut, no se sabe como varían en intensidad, ni como varían en función del 
tiempo. 

Cuando los bloques de la litósfera están sometidos a esfuerzos de compresión pueden dar 
lugar a desplazamientos correspondientes a fallas inversas, en cambio, si están sometidos a 
estados de tensión pueden dar lugar a desplazamientos de fallas normales y se crean los 
“graben”. 

Si miramos el perfil gravimétrico transcontinental, figura 2, (Introcaso? ), vemos que se 
puede dividir para simplificar su interpretación, en tres partes: 1) la que está al Oeste de la 
Ciudad de Chepes. 2) la que está comprendida entre Chepes y el rio Paraná, y 3) la que vá 
del río Paraná al Océano Atlántico. 

La primer parte presenta el aspecto típico de un borde continental, semejante al de la 
costa Sur de las Islas Sumatra y Java, que estudió y descubrió con todo detalle el Prof. 
Vening Meinesz.(*9) 

En nuestra Cordillera (a la latitud del perfil) la anomalía de Aire Libre alcanza el valor 
negativo de unos 100 miligales. Luego hacia el Este su valor medio se vá reduciendo hasta 
hacerse positivo en correspondencia con la Sierra de Chepes. 

La fuerte anomalía negativa indicaría que (si las fuerzas predominantes que actúan 
fuesen las que tratan de restablecer el equilibrio isostático), la Cordillera debería estar, en la 
actualidad, en un proceso de levantamiento. 

En el tramo medio del perfil, el detalle más importante es el valor positivo de la 
anomalía de Aire Libre en correspondencia con las Sierras de Córdoba. Si la anomalía es 
positiva, ¿a qué se debe el levantamiento de las Sierras? 

Sea como sea, tanto el detalle morfológico, como el dato geofísico son índices de 
inestabilidad tectónica, y por consiguiente de actividad sísmica. 
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El tercer tramo es el más uniforme. Relevamientos gravimétricos con mayores detalles, 


con estaciones más cercanas una de otra, podrían poner en evidencia probables líneas de 
fracturas, varias de las cuales son reconocidzs estratigráficamente. 


En resumen, los datos geofísicos aportan una información importante que corresponde 


tener en cuenta. Se hará cada vez más valioso a medida que se tengan más datos y se puedan 
correlacionar con otros de diferente origen. 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 
(12) 


19 Oct 1845 — 


5 Jun 1888 ' 


14 Dic 1924 
24 Ene 1927 
10 Nov 1932 
20 Jun 1933 
11 Jun 1934 


15 Jun 1934 


9 Ene 1935 

21 Ene 1948 
31 Jun 1960 
28 Oct 1965 


15 Oct 1968 


TABLA | 


Lista de los sismos representados en la figura 1 


16 - 


Lobos, Navarro, Chivilcoy, 
Costa del río Salado 


00 20 La Plata, Buenos Aires, 


2111 


S 30 


1353 


906 


307 


Campana, etc. 


ILAGAGSS 
IOSAAGSO 
309 619 
SLOP6SO 
33,3% 64,40 
Arias (Cba.) 


Tostado (Prov. de Sta. Fe) 


Goya-Esquina (Prov. 
Corrientes) 


Choele-Choel (Río Negro) 
Santa Rosa (La Pampa) 


Corzuela, Las Breñas 
(Chaco) 
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RIESGO SISMICO EN EL TERRITORIO ARGENTINO 


Fig. 1.-— Representación gráfica de algunos de los sismos más importantes que se produjeron en los 
últimos tiempos en zonas, del territorio argentino, consideradas como no sísmicas. 
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PERFIL TRANSCONTINENTAL 32S 


ALTIMETRIA 


ANOMALIA AIRE LIBRE 


_ TRAZADO LINEA 


Fig. 2. — Anomalía de Aire Libre a lo largo de un perfil transversal a la Cordillera de Los Andes a la 
latitud 329S que va del océano Pacífico al Atlántico (Introcaso 19723. 
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OBSERVACIONES SOBRE UN TEOREMA DE BRILLOUIN 


Por JORGE A. GUALA* 


RESUMEN 


Se examina el dominio de validez del teorema de L. BRILLOUIN sobre la distribución de la energía 
potencial entre dos partículas cargadas. Se establece que las consideraciones de simetría entre ambas cargas 
se cumplen en forma completa sólo si las dos partículas tienen igual masa. Finalmente, se demuestra que 
las partículas de masas My y M2 contribuyen a la energía potencial total con las cantidades (Q Q”/ry). 
(my/M;¡) y (Q Q/Ro) . (my/ma). respectivamente. 


SUMMARY 


The general validity of the BRILLOUIN'”s theorem about the potential energy distribution between 
two charged particles is considered. The symmetry considerations remain valid only if both masses are 
equal. We find the contributions of the particles of masses M, and M, to the total potential energy. They 
are (Q Q'/ro) . (my/M¡) and (Q Q'/ro) . (my/M)) respectively. 


En 1964 L. BRILLOUIN publicó un artículo! (L'énergie potentielle et samasse) en el 
que expone el problema de la distribución de la energía potencial entre dos cuerpos 
electrizados en interacción. Conocida la manera en que se reparte la energía potencial puede 
calcularse la distribución de la masa inercial ligada a esa energía. 

El razonamiento de Brillouin se funda en consideraciones electrostáticas, pues comienza 
por calcular la densidad de energía, u, correspondiente al campo de Coulumb, F, producido 
por un par de cargas, Q, Q”, aisladas en el espacio, en reposo relativo, a la distancia ry una 
de la otra. 

Así obtiene, para dicha densidad en un punto M, distante r de Q y r' de- Q”: 


02 2 ,2 > 
El [ qe E (1) 


8 7 87 E p? 1? p? 


La) 


S 7 
0 es el ángulo entre los vectores r y r. 


Seguidamente agrega: “Le premier terme (en Q?) intégré ¿ans tout l'espace, contribue a 
la masse propre M de la charge Q. Le second terme (en Q”) donne une contribution a la 


masse propre M” de Q”... Ler terme mixte (en Q Q”) S'integre aisément, reste toujours fini, 
et donne l'énergie potentielle 


* Centro Atómico. Bariloche. 


46 ANALES DE LA SOCIEDAD CIENTIFICA ARGENTINA 


Eso O: (Q) 


Cette énergie potentielle apparaít ainsi comme répartie dans tout l'espace, entre les charges 
Q Q” et tout autour. Le densité d'énergie (terme Q Q”) devient infinie sur Q et Q” (r =0 ou 
r” = 0). Il y a symétrie complete entre les deux charges et par conséquent il semble logique 
d'admettre que la masse de l'énergie potentielle se partage également entre les charges Q et Q”.” 

El último parágrafo del razonamiento de Brillouin merece ser examinado con deteni- 
miento, pues si bien existe simetría entre ambas cargas, ésta no es completa. En efecto, la 
aludida simetría es completa sólo desde el punto de vista electrostático, es decir, consideran- 
do el sistema ya conformado. Pero ocurre que para que las cargas Q Q” estén en reposo y 
separadas por la distancia ry fue necesario anteriormente trasladarlas a dicha configuración. 
Para poder efectuar este transporte de cargas se requiere la participación de un agente 
externo que provea una fuerza sobre cada partícula, que en todo momento sea igual y 
opuesta a la que ya actúa en virtud de la interacción electrostática. Es decir, hay que aplicar 
un par de fuerzas iguales y Opuestas a las culombianas, una sobre cada partícula, de modo 
que éstas puedan desplazarse con aceleración nula. Para cambiar la configuración del sistema 
de cargas, el agente externo debe realizar trabajo sobre cada partícula. Si se parte de una 
configuración inicial en la que ambas están en reposo 1elativo e infinitamente separadas, y se 
llevan a un estado final en el que también estén en reposo relativo pero distantes ry una de 
otra, entonces todo el trabajo realizado por las fuerzas exteriores se habrá transformado en 
energía potencial. 

Por ser la resultante de fuerzas cero, el centro de masa del sistema permanecerá en 
reposo, lo que obliga a que los desplazamientos de cada carga sean inversamente proporcio- 
nales a sus respectivas masas propias. Esto implica que las fuerzas aplicadas realizan trabajos 
diferentes sobre cada partícula Esto evidentemente rompe la simetría entre ambas cargas, 
excepto en el caso particular de ser sus masas idénticas. 

Como el trabajo efectuado sobre cada carga es inversamente proporcional a su respectiva 
masa propia, resulta que la energía potencial que adquiere cada una de ellas está también 
inversa de su masa, resultado enunciado por J. Palacios? (1968). 

La aparente simetría entre ambas cargas proviene de la forma de las ecuaciones 1 y 2. 
De esta última, por figurar sólo las cargas y la separación entre las mismas, podría sacarse la 
apresurada conlcusión de que la energía potencial de la carga Q” en el campo Q/ro 
engendrado por Q vale% Q'(Q/ro). La energía potencial de la carga Q en el campo de Q” 
sería, por idéntica razón % Q(Q”/ro). La energía potencial total valdrá la suma de ambas 
contribuciones, es decir Q Q”/ry y resulta evidente que ambas partículas contribuyen por 
igual. 

Es fácil mostrar que el anterior razonamiento no es correcto. Para ello expresamos la ec. 
1 en función de las distancias r,, r, de las cargas Q y Q” al centro de masa, respectivamente. 
Como 1, + r2 = To y My r, = M, t,, resulta: 


Eno = OA MO MS (3) 


my = M, M>/ (M, + M,) es la masa reducida. 
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Ahora expresamos la energía potencial total como suma de dos contribuciones, E; (r, ) y 
E>(r,). De acuerdo con la ley de Coulomb, E, y E, deberán variar como 1/r, y como 1/r, 
respectivamente. Como además se debe verificar la relación inversa entre energía potencial y 
masa, necesariamente E, deberá ser proporcionada a M>,/M,, en tanto que E, lo será a 
M/M). 

Con estas consideraciones podremos escribir: 


Eo LO POUM, PM.) € 3) = A O 


T r2 


donde C es una constante adimensional cuyo valor queda determinado por la condición E, + 
E, => Epot- 
En efecto, E, + E, como función de r, se expresa como 


C Eos (M, + M,) ) / My 


rt] 


Por comparación de esta -última ecuación con la ec. 3 resulta: 


C=M,/(M,+M,) 


con lo que 


ER 0 O a] 


E) = ——_—————- (M,/M, ) [m, / M, + M>)] 


Expresiones que, recordando que r, = r M2((M, + M)>) y r, = r M,/(M, + M)), toman la 
forma: 


Ae POMO EI DO MA) 


to To 


Las ecuaciones 4 y 4” muestran que si M, = M,, es decir m+ = 4M,,E,=E>»=%0 
Q'/To. 

For el contramor siiM. > Mom, = Mo? con lo que Ey =0 y E, = 0 Qro. 

Es decir, prácticamente toda la energía potencial está localizada en la más liviana de las 
partículas. 

Una comprobación de este resultado se obtuvo? en la evaluación teórica del potencial de 
ionización del Lj?. 
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